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A uniform asymptotic expansion of the time-harmonic antenna current with large parameter 
In [2kz/(ka)?] is given. As usual, the expansion is based on Hallén’s filament integral and is restricted 
by the condition z >> a but is otherwise uniform in kz. The use of averaging leads to the vanishing 
of the second term of the asymptotic series. Two different averaging procedures are used to derive 
two accurate formulas: The first one, which is very simple, is useful for practical applications. The 
second formula is of theoretical interest. An improved asymptotic series for large kz is derived to 
relate the averaging results to existing asymptotic results. Comparison of these two formulas and the 
exact antenna integral demonstrates that these formulas are more accurate than any other existing 
formulas for infinite cylindrical antenna current. The Hallén filament integral is found to approximate 
closely the exact integral even when the condition z >> a is violated. Discrepancies in the literature 


are discussed. 


1. INTRODUCTION 


It is interesting to note that early authors hesi- 
tated to consider an infinite antenna. Hallén [1948] 
had to justify its consideration by saying, ‘“‘An 
infinite antenna does not exist, but the behavior of 
the outgoing traveling current waves on an antenna 
of finite length must be the same as if it were 
infinite ....°’ The need for good infinite antenna 
solutions is universal: Hallén [1962] and King [1956] 
used an infinite antenna solution to illustrate the 
difference between the antenna and the transmis- 
sion line solutions. In antenna handbooks [Rispin 
and Chang, 1988; King and Smith, 1981; Vain- 
Shtein, 1969] the finite antenna solution is con- 
structed by means of an infinite antenna solution. 

The earliest treatment of an infinite cylindrical 
antenna was given by Stratton and Chu [1941]. 
Other important works on antenna admittance and 
feedgap problems were given by Brillouin [1943], 
Infeld [1947], Albert and Synge [1948], Synge 
[1948], Wu and King [1959], Duncan [1962], Chen 
and Keller [1962], Fante [1966], and Latham and 
Lee [1969]. Two Russian articles treated an infinite 
length antenna [Vainshtein, 1961b] and a finite 
length antenna [Vainshtein and Fock, 1968]. The 
first correct asymptotic solution valid to the leading 
(zero) order was obtained by Vladimirskii [1944], 
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and to the first order a little later by Hallén [1948]. 
Numerous other authors [Vainshtein, 1959, 1961; 
Kuehl, 1960; Wu, 1961; Kunz, 1963] have discussed 
asymptotic expansions of antenna current. 

Papas [1948] obtained an asymptotic expansion 
of the magnetic field away from the antenna for an 
infinite cylindrical antenna and correctly conjec- 
tured that the form of the solution for the antenna 
current would be similar. Northover [1958] dis- 
cussed asymptotic solutions for the magnetic field 
and the current of an infinite antenna, and pointed 
out a difference of a factor of 2 when comparing the 
magnetic field away from the antenna and the 
antenna current. Latham and Lee [1970] obtained 
the time domain magnetic field and showed graph- 
ically a factor of 2 transition from the antenna 
current to the magnetic field away from the an- 
tenna. 

The first uniform formula for infinite antenna 
current was obtained in the work by Vainshtein 
[1961] in the Russian literature. Shen et al. [1968] 
also gave an interpolation formula, which will be 
referred to as the SWK interpolating formula. 
Weinstein’s formula has been used in the antenna 
literature for constructing the scattering antenna 
current excited by a grazing incident plane wave 
[Chang et al., 1978]. 

Averaging has been shown to be a powerful way 
to construct uniform formulas [Chen, 1983, 1987; 
Chen and Warne, 1991]. The gist of averaging has 
been discussed in the treatment of time domain 
antenna current [Chen and Warne, 1991]. Averag- 
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ing is used here to derive two formulas for time- 
harmonic current: One formula is very simple; the 
other shows the connection to Weinstein’s formula 
and the more general solution for loaded antennas 
[Chen and Warne, 1992]. These uniform formulas 
are not only the most accurate formulas, but they 
also can be used to understand the small differences 
between the antenna solution (in the limit of van- 
ishing ka) for moderate length and transmission line 
solutions. Pocklington’s [1897] sinusoidal distribu- 
tion is, of course, basically a transmission line 
solution. The great difference between the antenna 
and transmission line solutions only applies to long 
antennas. The time-harmonic and time domain be- 
havior of a long antenna in a dissipative medium, 
however, also resembles a transmission line solu- 
tion [Chen, 1987]. 

Finally, some explanations on the asymptotic 
solution of infinite antenna current are helpful in 
understanding the literature, since the infinite an- 
tenna is such a basic problem in electromagnetic 
theory. The main objective, however, is to derive 
accurate formulas for the time-harmonic antenna 
current and show the accuracy of these new formu- 
las. 


2. DERIVATION OF A SIMPLE FORMULA 
FOR THE TIME-HARMONIC ANTENNA CURRENT 


The time-harmonic current on an infinite tubular 
antenna (with radius a) driven at a delta gap by a 
voltage — V is (with e~' time dependence) [Hallén, 
1948, 1962; Wu, 1969] 


2ikV 
I(z) =e 
no 


eifz dé 
j. (k? — €?)ariTolatk? — €7)" Lg Lat? — €7) 12)" 


(1) 


where 19 = ({9/€9) 12 is the intrinsic impedance of 
the medium and k = o/c. The usual contour of 
integration C is indented below the branch point at 
k but indented above the branch point at —k (Figure 
1). When the contour is so chosen J(z) in (1) gives 
the sum of the exterior and interior currents. Be- 
cause the current (1) is even with respect to z, we 
can confine our attention to z > 0. The interior 
current arises from the residue contribution of the 
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Fig. 1. 


The contour of integration is indented below the 
branch point at k but indented above the branch point at —k. 
When the contour is deformed to the upper half plane, the 
interior current arises from the residue contribution of the zeros 
of Jola(k? — €2)"2] lying on the positive imaginary axis (shown 
as ‘‘x’’), When the exterior antenna current is the subject of 
study, we do not include these residue contributions. 


poles (zeros of Jo) lying on the positive imaginary 
axis. When the exterior antenna current is of pri- 
mary concern as in this paper, we do not include the 
residue contributions of all poles on the imaginary 
axis [Hallén, 1962]. 

To obtain an asymptotic solution for large z, it is 
permissible to approximate the Bessel functions by 
the leading terms of their small-argument expan- 
sions [Vladimirskii, 1944]. However, it should be 
emphasized that Hallén gave great significance to 
such an expanded form by interpretating it as a 
simple physical space kernel. The solution obtained 
can be used for wires of general cross section and is 
also free from the logarithmic singularity associated 
with the gap capacitance. Any realistic feedgap can 
be defined approximately as the difference between 
this idealized configuration and the actual configu- 
ration, at least for thin antennas. An example of 
such a treatment can be found in the discussion of 
effective impedance of a bolt load on a narrow 
aperture [Warne and Chen, 1992]. Another ideal- 
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ized configuration was treated by Latham and Lee these terms in the denominator of (6) are small 
[1969]. anyway. Thus (6) becomes 
After the incorporation of the small-argument 


expansions, (1) is reduced to [Hallén, 1948; _.f? Jeep i/2 
Vladimirskii, 1944] G(Bp, kz) =i { dne™" > = ro 
1+i- 
2ikV 2 
I(z) ~ — 
No 1 
ig ; 3 
-{ ——o+——_ pp -Ing+in(1+i2) 21% 
A : 2 2 
C (2 - e-2 In (5) ~ In (k? ~ £2) -2y + a 
1 
; ey fae Vee (7) 
2) pp-inn in (i e2) 1 


Equation (2) can be written as 
A subsequent integration by parts leads to 


VK 
I(2) a9 e™G(Bp, kz), (3) 
0 ‘00 
G(Bp, kz) | dne“™ 
where 0 
eixkz dy n 7 
(Bp, ke) = -1 | 7 ; Bp-In(n) +In [1 +45] -i > 
oA(14d Bp —In (ied *lo |} ——_—_—_———————< |, (8) 
2 2 7 39 
(4) Bp —1In(m) +1n Pa are 
with 


Two distinct averaging procedures will be used 
&=K(1+ x), Bp = —2 In (ka) — 2y + In 2. (5) below. First, a simple averaging of In (») over the 


weight e~ ™ is defined as 
Because of the small-argument expansion of the 


Bessel functions, extraneous poles located at x = 20 
—1 + [2e% + 112 have been introduced [Hallén, i) 
1962; Wu, 1969]. Here we do not go into the 


dn In (n)e~™ 
2 = —In (kz) — y. (9) 


discussion of this matter but simply note that the co ae 
contour C is to be above these poles for z > 0. Thus { dne 
G(Bp, kz) : ; F ‘ 
Using (9) in (8) and the exponential of (9) for 7 in 
iF dne~™® 1 In (1 + inf/2), we obtain 
Sj es 
0 a .1\| 3m 
a(re2) Bp fal i2)] oe Bp + y—In2 + In (2kz + ie-%) — 1 
G(Bp, kz) ~ § |] AAA a 
1 37 
= ¥ (6) Bp + y —In2 + In(2kz + ie”) + i — 
—In 1+i a Vhs i me 
Bp n 2 2 
WT 
Equation (6) can be simplified if we change the signs —21n (ka) — y + In(2kz + ie~”) —i 7 
of the denominator terms In (1 + in/2). The approx- = i In| ————_—_——————_————__|__ (10) 


imation in general requires |B,| >> 1 (Note further- 


WT 
—21n (ka) — y + In(2kz + ie~%) + i — 
more when kz is large In (1 + i7/2) > 0.) so that ple Sia ea ae ks 2 


652 


3. DERIVATION OF AN ALTERNATIVE FORMULA 


An alternative averaging procedure can be de- 
fined as 


ie dye" [i (In (: + i3)| 
; 2 
[ dne 
0 


= -In (kz) — y — e!?®E, (—i2kz). (11) 


Here E, is an exponential integral. Use of (11) in (8) 
gives 


Gap, te) ~ in| 12 
ener = erage" ue 
where 

2kz Me Ss os : 
Q(z) = In (kay? —-yti 2 +e ‘E,(—i2kz). (13) 


Note that E, can be expressed as sine and cosine 
integrals: 


E\(—i2kz) = —Ci(2kz) — isi(2kz), (14) 


where si(2kz) = Si(2kz) — 2. 


4. SUMMARY OF EXISTING FORMULAS 
FOR ANTENNA CURRENT 


Frequently used formulas in the literature are 
compared to the new formulas (10) and (12). In the 
next section they are compared to the numerical 
integration values of the exact solution (1) and the 
Hallén filament approximation (6). 

The Hallén first-order formula (good for two 
terms) is the earliest first-order asymptotic result 
for large parameter In [2kz/ (ka)”]. Let © be given 
by (O is asymptotically equivalent to Q(z) as kz > 


oo.) 


a 2kz ae 7 
=In (ka) yti 5" 
Then, Hallén’s formula (with G(8,, kz) defined in 
(3)) is [Hallén, 1948] 


a (15) 


2a 
G(Bp, kz)~—- (16) 
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Equation (16) is also discussed by Vainshtein 
[1959]. 

Two approximate formulas for infinite cylindrical 
antennas are described by Shen et al. [1968]. A 
truncation method is used to derive an approximate 
formula for large kz: 

Q-in 
Q+ in| 


This formula is also derived in the limit as kz > ~ by 
averaging [Chen, 1983] and will be shown later to be 
a first-order asymptotic solution for a different 
choice of a large parameter. Equation (17) will be 
referred to as the SWK (Shen, Wu, and King) 
formula. The second formula by Shen et al. [1968] is 
obtained by interpolating (17) and the antenna input 
admittance derived by Wu [1961]: 


G(Bp, kz) ~iln (17) 


G(Bp, kz) ~i 


—2In (ka) — y + In{kz + [(kz)? + e277] '} - is 

=| ——_————_—_—_—_—_———cxK«“ |, 
3a 
—2In (ka) — y + In {kz + [(kz)? + e727]! + iF 


(18) 


Equation (18) is the SWK interpolating formula. 
Note that this formula has been used for construct- 
ing finite antenna solutions [Rispin and Chang, 
1988; King and Smith, 1981]. 

Finally, (12) with Q(z) > 1 gives 


2a 
G(B,, kz) eee pax: 


yon a 


Equation (19) was previously derived by Vainshtein 
[1961a] using a variational technique and has been 
used widely for constructing finite antenna solu- 
tions [Rispin and Chang, 1988]. This result could 
have been obtained by applying averaging directly 
on (4). (First, In (1 + y/2) is neglected, and then 
In () is replaced by its average defined over the 
loop integral with {e2/[(1 + x/2)]} as the weight.) 
When kz — ~, (17) can be shown to be 


ar 

307)" 
Chen and Wu [1969] observed that the expansion of 
an identical formula is accurate to the leading order 


2a 
G(Bp, ke)~ [I (20) 
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in 1/0, but actually, the leading term of (20) is good 
to the first order. Equations (10), (12), (17), and (19) 
reduce to (16) as kz > ™. 

When kz — 0, (10) and (12) with (3) become 


10) i | —In (ka) — y 
n 


=In (ka) — y + im el} 


Vso 


and (19) approaches 


10) ow 1 
me i (ka) — y + Fi 


Note that when ka is very small (—In (ka) > 1), (21) 
can be expanded to yield (22); when ka is not very 
small, (21) gives more accurate numerical values. 
Equation (22) is the antenna input admittance de- 
rived by Tai [1955], Hallén [1962], Vainshtein 
[1959] and Borgnis and Papas [1955]. Note that as 
ka — 0, (18) becomes 


(22) 


7 oT 
—In (ka) - y+i-—-it= 


10) i 4 2 
—~— in rrr aie 
am a ye ee ee 
n(ka)—ytit+is 

7 1 


—. (23) 
we [-in (ka) - iG] 


It is not surprising that (23) differs from the correct 
result (22) because only the conductance is used in 
the interpolation. 


5. NUMERICAL RESULTS 


Numerical integration of G, defined by (3), for the 
exact solution (1) and the Hallén filament approxi- 
mation (6) is facilitated by subtracting the left-hand 
side of (24) below and then adding its explicit form 
(the right-hand side of (24)). This facilitating slowly 
convergent integral is 


(kz)! d 
Go(Bp, kz) =i { : 


0 n 
1+i- 
1 


3a 


” 
-Ingn+infl+i-]|+i— 
Bp nN ( 2) Ls 


3 
A-- 
1 2 
— | =i In : 
Iny +In ee i Pa 
—In i-|-it- 7 
Bp Inn 2 2 2 
(24) 
where 
1 1 
= © arian (7) +1 
7 2kz 
2/1 
-i- zi (1 + 4k2z?) — In (ka) — y|. (25) 
7 


Table 1 shows a comparison of the exact solution 
(1) and (3), the Hallén filament approximation (6), 
the simple averaging formula (10), the alternative 
averaging formula (12), the SWK formula (17), the 
Hallén first-order formula (16), the Weinstein for- 
mula (19), and the SWK interpolating formula (18) 
for various ka and kz values. Numerical values are 
determined to greater accuracies than those shown 
in the table. Of particular interest is that the Hallén 
filament approximation is in agreement with the 
exact solution; even when ka = kz = 0.1, the real 
part of G has a 0.3% error and the imaginary part of 
G has an 8% error. This accuracy is consistent with 
the time domain case discussed by Chen et al. 
[1992]. The simple averaging, the alternative aver- 
aging, and the Weinstein formula are found to be 
comparable in accuracy. The SWK interpolating 
formula has somewhat larger errors because, as 
shown in (23), the imaginary part of the antenna 
input admittance is not used in the derivation of the 
formula. SWK and the Hallén first-order formulas 
are derived for kz > ©; however, remarkably, when 
kz = 0.5, both the real and imaginary parts are quite 
accurate. As expected, the SWK formula is slightly 
more accurate than the Hallén first-order formula, 
especially when ka is not very small. 

Figures 2 a and 2b show, for ka = 0.5 and ka/10 = 
kz = 1000, the real and imaginary parts of G, 
respectively, for all eight entries discussed in Table 
1. Note that even at this large ka value (for realistic 
antennas) the Hallén filament approximation has 
only 1% error in the real part when kz = ka/10; the 
imaginary part given by the dotted curve, however, 
has a large (55%) error. A realistic antenna feedgap 
is probably of the order of ka in width; in such a 
case the discrepancies of the current obtained from 
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TABLE 1. Comparison of Numerical Integration of the Exact Solution (1) and (3) and the Hallén Filament Approximation (6), 
Weinstein Formula (19), and the SWK 
Exact Solution Hallén Filament Approximation Simple Averaging 
kz Real Imaginary Real Imaginary Real Imaginary 
0.001 (=ka) 
0.001 0.47952347 —0.12612700 0.47952406 —0. 12549869 0.4607 —0.1100 
0.010 0.47848421 —0.11681319 0.47848441 —0.11679040 0.4612 —0.1089 
0.100 0.47045314 —0.09447538 0.47045318 —0.09447468 0.4638 —0.0983 
0.500 0.44979931 —0.06892960 0.44979932 —0.06892952 0.4514 —0.0679 
1.000 0.43544673 —0.05833887 0.43544674 —0.05833884 0.4354 —0.0557 
2.000 0.41905322 —0.04972731 0.41905322 —0.04972730 0.4176 —0.0474 
5.000 0.39679600 —0.04164043 0.39679600 —0.04164043 0.3948 —0.0403 
10.000 0.38075004 —0.03733794 0.38075004 —0.03733794 0.3788 —0.0364 
0.010 (=ka) 
0.010 0.72067080 —0.28649826 0.72072689 —0.28021016 0.6640 —0.2354 
0.100 0.70268351 —0.21351438 0.70269988 —0.21328330 0.6740 —0.2137 
0.500 0.65901916 —0.14721664 0.65902299 —0.14719598 0.6574 —0.1446 
1.000 0.6297 1087 —0.12134681 0.62971277 —0.12133866 0.6269 —0.1153 
2.000 0.59700466 —0.10055131 0.59700555 —0.10054792 0.5923 —0.0951 
5.000 0.55383346 —0.08099890 0.55383377 —0.08099778 0.5488 —0.0775 
10.000 0.52360703 —0.07057560 0.52360716 —0.07057510 0.5190 —0.0681 
0.100 (=ka) 
0.100 1.31231430 —0.75973227 1.31667978 —0.69527401 1.1343 —0.6711 
0.500 1.17988440 —0.46218380 1.18179344 —0.45447198 1.1468 —0.4412 
1.000 1.09721339 —0.36082412 1.09816573 —0.35805664 1.0748 —0.3346 
2.000 1.00902834 —0.28270580 1.00945101 —0.28167250 0.9878 —0.2603 
5.000 0.89858823 —0.21132588 0.89871995 —0.21102175 0.8803 —0.1968 
10.000 0.82528611 —0.17439130 0.82533908 —0.17426576 0.8099 —0.1641 


the Hallén filament approximation have a 10% error 
in the imaginary part. We can thus assume that the 
Hallén filament approximation typically has small 
error when applied to realistic antennas. The two 
averaging formulas have small errors, and the sim- 
ple averaging is closest to the exact solution. The 
SWK interpolating formula gives a fairly accurate 
real part for this case. The Weinstein formula does 
not give very good numbers for this large ka. 
Figures 2c and 2d show, for ka = 0.1 and ka/10 = 
kz = 1000, the real and imaginary parts of G, 
respectively, for all eight entries. For the useful 
range of kz = ka the two averaging formulas are 
the closest to the exact solution and the Weinstein 
formula is the next closest. The SWK interpo- 
lating formula gives fairly accurate numbers for 
the real part. Figures 3a and 3b show, for ka = 0.01 
and ka/10 = kz = 1000, the real and imaginary 
parts, respectively, for all eight entries. The two 
averaging formulas and the Weinstein formu- 
la are comparable in accuracy. The averaging 
formulas are more accurate than the Weinstein 


formula for kz somewhere between 0.1 and 1, 
and vice versa for kz not in this range. Figures 3c 
and 3d show, for ka = 0.001 and ka/10 = kz 
= 1000, the real and imaginary parts, respective- 
ly, for all eight entries. The same comment that 
applied to Figures 3a and 3b also applies to 
Figures 3c and 3d. 

In conclusion, it is interesting to note that the 
simple averaging formula (10) is the simplest and 
also gives the most accurate approximation to the 
exact solution throughout the parameter range stud- 
ied. The alternative averaging formula (12) gives 
almost as accurate numbers as the simple averag- 
ing. Next in order of accuracy is the Weinstein 
formula (19). The SWK interpolating formula (18) is 
somewhat less accurate, giving too large an imagi- 
nary part for G (the imaginary part of G is negative). 
For asymptotic formulas, the SWK is slightly more 
accurate than the Hallén first order, but both have 
very small errors when kz = 0.5. Finally, the Hallén 
filament approximation has essentially no error for 
practical antenna applications. 
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the Simple Averaging Formula (10), The Alternative Averaging Formula (12), SWK (17), the Hallén First-Order (16), the 
Interpolating Formula (18) 


Alternative Averaging SWK Hallén First Order Weinstein SWK Interpolating 


Real Imaginary Real Imaginary Real Imaginary Real Imaginary Real Imaginary 
0.4608 —0.1096 0.8110 —0.1614 0.8520 —0.1905 0.4676 —0.1155 0.4600 —0.0526 
0.4610 —0.1069 0.6347 —0.0997 0.6551 —0.1103 0.4680 —0.1126 0.4595 —0.0525 
0.4586 —0.0926 0.5193 —0.0670 0.5306 —0.0717 0.4658 —0.0976 0.4544 —0.0514 
0.4440 —0.0693 0.4601 —0.0527 0.4680 —0.0555 0.4509 —0.0728 0.4352 —0.0472 
0.4314 —0.0585 0.4385 —0.0478 0.4453 —0.0502 0.4378 —0.0613 0.4199 —0.0439 
0.4160 —0.0495 0.4187 —0.0437 0.4247 —0.0456 0.4218 —0.0517 0.4031 —0.0405 
0.3944 —0.0412 0.3951 —0.0389 0.4002 —0.0405 0.3995 —0.0429 0.3815 —0.0363 
0.3787 —0.0369 0.3790 —0.0358 0.3834 —0.0371 —0.3832 —0.0382 —0.3664 —0.0335 
0.6646 —0.2309 1.1057 —0.2932 1.1982 —0.3987 0.6800 —0.2569 0.6784 —0.1137 
0.6650 —0.1996 0.8110 —0.1614 0.8520 —0.1905 0.6830 —0.2224 0.6676 —0.1102 
0.6415 —0.1456 0.6796 -0.1141 0.7045 —0.1282 0.6606 —0.1613 0.6282 —0.0977 
0.6180 —0.1202 0.6347 —0.0997 0.6551 —0.1103 0.6359 —0.1323 0.5976 —0.0885 
0.5887 —0.0991 0.5951 —0.0878 0.6121 —0.0960 0.6047 —0.1081 0.5647 —0.0791 
0.5478 —0.0793 0.5495 —0.0750 0.5629 —0.0809 0.5610 —0.0855 0.5240 — 0.0682 
0.5186 —0.0690 0.5193 —0.0670 0.5306 —0.0717 0.5299 —0.0738 0.4964 —0.0612 
1.1216 -0.6171 1.6726 —0.6080 1.8272 — 1.1867 1.1048 —0.8058 1,1997 —0.3414 
1.0981 —0.4329 1.2354 —0.3604 1.3540 —0.5280 1.1485 —0.5760 1.0875 —0.2841 
1.0464 —0.3436 1.1057 —0.2932 1.1982 —0.3987 1.1095 —0.4492 1.0050 —0.2445 
0.9760 —0.2696 0.9986 —0.2415 1.0704 —0.3105 1.0374 —0.3417 0.9205 —0.2065 
0.8773 -0.2011 0.8831 —0.1906 0.9350 —0.2320 0.9268 —0.2440 0.8220 —0.1657 
0.8088 —0.1661 0.8110 —0.1614 0.8520 —0.1905 0.8491 —0.1958 0.7587 —0.1417 

ka=0.5 
3.0 —— Exact 
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a Alternative Averaging 
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1) 
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Ps 
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Fig. 2a. The real part of the antenna current normalized to give G = Io/Ve, which is (3). 
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Fig. 2b. The imaginary part of the antenna current normalized by use of (3). 


6. AN ASYMPTOTIC EXPANSION OF ANTENNA 
CURRENT 


In this section we derive an asymptotic expansion 
of antenna current and relate our expansion to 
existing works. Although the finite antenna solu- 


Fig. 2c. 


tion had previously been given by Hallén, let 
us comment that for infinite antennas the first 
zero-order formula for input impedance and the 
asymptotic formula, for kz — ~, together with 
proper interpretations of the meaning of these 
formulas, were apparently first given by 
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Fig. 2d. The imaginary part of the antenna current normalized by use of (3). 


Vladimirskii [1944]. He also briefly discussed 
the logarithmic singularity of the gap capaci- 
tance. 

In order to demonstrate the use of averaging 
for deriving an asymptotic series, we consider the 
limit as kz — ©, We can thus approximate 


0.4 


03 
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Fig. 3a. 


(1 + én/2) in (8) by 1 and expand the logarithmic 
factor in the integrand in a Taylor’s series in 
In (7m) at the expansion point x9. When the 
expansion point is chosen as x9 = —In (kz) — y 
given by (9), the second term of the series vanishes, 
and the resulting expansion for G(B,, kz) is 
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QO - in where 
G(B,, kz) ~ il 
(pst) In OQ+in 20 
Pn =kz { [In (1) + In (kz) + y]"e~ 2" dn. (27) 
oa na yn ||? The quantity p,, can also be determined by differ- 
Q+ a- 26 n 
n=2" (+ im" bn) i a entiating the following generating function: 


ka=0.001 


o 
D Seveeadee: Simple Averaging . 
~% = Altermative Averaging 

0.35 SWK. 

0.30 

0.25 

o.0oof 0.001 0017 o17 17 107 1007 1000 

kz 


Fig. 3c. The real part of the antenna current normalized by use of (3). 
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Fig. 3d. The imaginary part of the antenna current normalized by use of (3). 


n! 
n=0 


g(s) = kf (nkze")Se~" dn = eYT(1 +s) = Ba is”. 
0 
(28) 
The leading coefficients are 
po=l, 
pi =9, 
pz = £(2) = 77/6 = 1.644934, 
p3 = —2¢(3) ~ —2.404114, 
p4 = 3774/20 ~ 14.611364, 
ps = —24£(5) — 1077/3 £(3) ~ —64.432353, 
Po = 40[(3)]* + 6127 9/168 ~ 406.873470, 
p7 = —720£(7) — 842r7¢(5) — 2124/2 £(3) 
=~ —2815.131424, 
and 
pg = 2688L(3)¢(5) + 5607/3 [£(3)]? + 1261 78/720 
~ 22630.60734. 


The leading term in (26) was given in (17) and 
discussed before. 


Kunz [1963] gave an asymptotic expansion of 
antenna current in terms of the sum of two series 
with one in inverse powers of 9 and another like 
series divided by kz. He pointed out that the first 
series is dominating and computed the coefficients 
of a few leading terms of both series. He also 
claimed that the leading term of his first series is 
more accurate than Hallén’s first-order formula 
(16). It is clear that Kunz’s leading term has a 
smaller absolute value in the denominator than that 
of Hallén’s and thus would make the numerical 
value farther away from the numerical integration 
value of the exact integral. What Hallén did was to 
choose the large parameter such that the second 
term of the series vanishes, which is in the same 
spirit as the averaging method. 

Hurd [1966] identified the general coefficient as 
essentially the nth derivative of the Gamma func- 
tion evaluated at unity. Note that the transforma- 
tion from this general coefficient to the coefficient of 
Kunz’s series is not readily available. In the present 
asymptotic representation a similar general coeffi- 
cient is used directly for the expansion, and thus the 
extra transformation has been eliminated through 
the choice of a large parameter. 

Furthermore, the use of averaging eliminates the 
first-order term, and thus the leading term is good to 
the first order (two terms). The first term alone is a 
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reasonable approximation even if kz is of the order 
of unity. Of course, to apply this result as an 
approximation to the original exact integral (1), the 
restriction of the Hallén filament approximation 
z => a must still be observed. Numerical compari- 
sons of (1) and (2) or (6) in the last section has 
relaxed this condition considerably. 

Kuehl [1960] applied a saddle point integration to 
evaluate infinite antenna current for kz > ~. It is 
interesting to note that the nature of the integral (2) 
is such that no large parameter is available in the 
exponent, and thus one cannot approximate the 
integration along the steepest descent path by a 
Gaussian [De Bruijn, 1981], as was done by Kuehl 
[1960]. 


7. CONCLUSIONS 


Two averaging formulas are derived. The simpler 
‘*simple averaging formula’’ (10) is slightly more 
accurate than the alternative averaging formula (12) 
when compared to the numerical integration of the 
exact solution (1) and (3). This simple averaging 
formula is useful for all practical applications. The 
alternative averaging formula is of theoretical inter- 
est in general antenna solutions. The Hallén fila- 
ment approximation (6) is also compared to the 
exact solution and found to have typically small 
error when applied to realistic antennas. The Wein- 
stein formula (19) has a slightly larger error than the 
two averaging formulas for a wide range of param- 
eters. The SWK interpolating formula (18) has the 
largest error of all the uniform formulas but is still 
quite accurate. The two asymptotic formulas: the 
Hallén first order (16) and SWK (17) are accurate 
even when kz = 0.5; SWK is slightly more accurate 
than the Hallén first order. An asymptotic expan- 
sion of antenna current for kz — ~ is derived, and 
the leading term is shown to give rise to the SWK 
formula. 
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